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The EXOTIME project: 

a status report on PG 1325+101 (QQ Vir) 
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After the discovery of V391 Peg b, the first planet detected around a post Red Giant phase star (Silvotti 
et al. 2007) , the EXOTIME (EXOplanet search with the Timing MEthod) project is focused on the search 
for new planets with similar characteristics. The aim of the project is to organize a global observing network 
to collect as much data as possible for a sample of five subdwarf B (sdB) stars and share them in order to 
obtain a more precise analysis. These evolved pulsators may have extremely regular oscillation periods. This 
feature makes these stars suitable to search for planetary companions with the timing method as in the case 
of pulsars. In this contribution we present the project and some preliminary results for the star PG 1325+101 
(QQ Vir) after the first two years of activity. 
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1 Introduction 

The EXOTIME project (EXOplanet search with the 
Timing MEthod) is a coordinated observing program 
aimed at the search for substellar companions around 
pulsating subdwarf B stars and the derivation of evo- 
lutionary timescales. The required observations span 
over a very long time (of the order of years), so this 
method is quite expensive in terms of observing time, 
but with the advantage to be sensitive to wide orbits 
and relatively low masses (down to ^ AIj„p). 

Using the timing method we can measure both the 
variations of the pulsation periods and the phase varia- 
tions; the latter potentially allow us to detect the pres- 
ence of substellar companions as in the case of V391 
Peg b (Silvotti et al. 2007), the first planet discovered 
with this technique around a pulsating star. The dis- 
covery of V391 Peg b has raised the interest to inves- 
tigate evolved planetary systems beyond the main se- 
quence and beyond the red giant branch. The orbital 
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distance of this planet, lower than 2 AU, suggests that 
this planet may have "survived" to the RG expansion 
of the parent star. Recently, two substellar compan- 
ions have been detected orbiting the sdB eclipsing bi- 
nary HW Vir (Lee et al. 2009), suggesting that substel- 
lar objects might be a relatively common phenomenon 
around sdB stars. 

Further goals of EXOTIAIE are the characteriza- 
tion of the targets using asteroseismic methods and 
the measurement of the secular variation of the oscil- 
lation period (the so called P-dot, P), giving a further 
example of the synergy between asteroseismology and 
the search for exoplanets, already known in the case of 
solar-like stars. Measuring P allows a precise determi- 
nation of the evolutionary status of a star and can help 
the identification of the pulsation modes. 

As a final goal, EXOTIME wishes to improve our 
understanding of the formation and evolution of the 
sdB stars. The formation processes of sdB stars still 
represent an unclear topic in stellar evolution. Different 
scenarios have been proposed, both for single stars and 
binaries. The presence of a secondary body such as a 
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planet, in particular for single stars, has been suggested 
to play a role in this process, enhancing the mass loss 
near the RGB tip (Soker 1998). 

2 The Method 

As in the case of the pulsar timing, we use the os- 
cillation periods as a clock in order to detect all the 
possible variations in the travel time of the photons. 
These variations could be attributed to the secular vari- 
ation of the oscillation period, P, or to a wobble of 
the sdB location due to the presence of a perturbing 
body, such as a planet. These variations are easily de- 
tectable through the 0-C diagram, in which the theo- 
retical [Calculated) expectations of a particular quan- 
tity are compared with the Observed ones. Once one 
obtains a sufficient number of data, it is possible to 
compare the mean phase of a long monitoring period 
with the phases at different time steps. If changes oc- 
cur in the pulsation period, they can be detected and 
identified according to the distribution of the points 
in the graph. When a pulsation period changes lin- 
early in time, the 0-C diagram shows a parabolic shape 
caused by the evolutionary timescale. The presence of 
a companion is revealed by a sinusoidal trend which 
means cyclically advanced or delayed timings of the 
maxima (or minima), due to the motion around com- 
mon barycenter. 

The drawback of timing is the need of long-term 
monitoring, which imposes to collect data for many 
years. Anyway this technique is somewhat complemen- 
tary to the radial velocity and transit methods, since 
these methods can hardly provide information on stars 
with small radii and hardly detect planets in wide or- 
bits and relatively low masses (down to ~ Mj„p). 

3 Targets 

The targets of EXOTIME (see Table [J) were selected 
following different criteria. We excluded known bina- 
ries and stars with an IR excess from the 2MASS data. 
Other important features are stable amplitude spectra 
in terms of phases and amplitudes and simple spectra, 
ideally with only 2 or 3 frequencies. A low number of 
frequencies allows to resolve the frequency spectrum 
even in short runs while, on the same time, it gives the 
opportunity to obtain independent O-C plots from each 
individual frequency. Finally we gave higher piority to 
bright targets, high pulsation amplitudes and high ob- 
servability, considering that the available telescopes are 
mostly in the northern emisphere. 

4 Observations 

Since 2008 the EXOTIME collaboration collected a 
large amount of data from several observing sites cov- 



Table 1 Targets of the EXOTIME program. Ref- 
erences: (a) Ostensen et al. 2001b; (b) Dreizler et al. 
2002; (c) Kilkenny et al. 2006; (d) Silvotti et al. 2002; 
(e) Ostensen et al. 2001a. * Hybrid pulsators. 



Object name 


ref. 


B 


Main 


Amplitude 






mag 


Period [s] 


[mma] 


HS 0444+0458 


a 


15.2 


~137 


~11 


(V1636 Ori) 










HS 0702+6043 


b 


14.7 


~360 


~22 


(DW Lyn)* 










EC 09582-1137 


c 


15 


~136 
~151 


~8 
~7 


PG 1325+101 


d 


13.8 


~138 


~26 


(QQ Vir) 










HS 2201+2610 


e 


14.3 


~350 


~10 


(V391 Peg)* 











ering the longitude range between Eastern Europe to 
Western North America, equipped with telescopes hav- 
ing typical apertures of 1-2 meters (up to 4 m) . The ob- 
servations are performed in B filter and, according to 
the pulsation period of the target and the magnitude, 
the time sampling is set to obtain at least 5-6 points 
per cycle. Because of the different characteristics of the 
instruments, it is necessary to weight the incoming data 
in order to obtain a coherent data set. 

The list of the observations performed during the 
first two years of EXOTIME is available from a ded- 
icated web site {www. na. astro. it/'^ silvotti/ exotime/). 
The web site provides informations on the target stars, 
observing instructions and scheduled observations, use- 
ful to optimize the planning of the observations, and 
the coverage of the targets. 

5 Status of the project 

The EXOTIME program carries on with its activity of 
collecting data for each target star. Preliminary analy- 
sis of HS 0702+6043 and HS 0444+0458 are available 
in Schuh et al. 2010, while an updated phase analysis 
including the new data on V391 Peg will be performed 
in the next months. The data set collected up to now on 
EC 09582-1137 is still too poor for our goals. Here we 
present the preliminary analysis of PG 1325+101 (QQ 
Vir) using part of the available data. Since we know 
that many of the sdB stars show amplitude variations 
(and sometimes phase variations) on various time scales 
from years to days, a preliminary check on the stability 
of the oscillations is required. 

5.1 PG 1325+101: preliminary analysis 

As a first step we can exclude the presence of a cooler 
stellar companion, studying the Spectral Energy Dis- 
tribution. Figure [T] shows the comparison between a 
black body (which provides a good approximation of 
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Fig. 1 Spectral Energy Distribution of PG 
1325+101. The filled red circles are 2MASS JHK 
and BVR magnitudes from the literature. There is no 
evidence for a cool stellar companion to the sdB. 



the energy distribution of a sdB star) and the 2MASS 
JHK and BVR magnitudes. The plot doesn't show any 
infrared excess. 

After this check we can start to analyse the incom- 
ing data set for this target. Generally the data are re- 
duced directly by the observers, but since the outputs 
can be expressed in different ways, it is very important 
to turn them into a standard format. The time series 
are then barycentrically corrected and the differential 
photometric flux is expressed in mmi (milli-modulation 
intensity units, 1 mmi=0.1%). Then we calculate the 
uncertainty for each data point (following the approach 
of Silvotti et al. 2006) in order to fix the weights to the 
different data sets. The available data are shown in Ta- 
bleland they sample the seasons 2008 and 2009 from 
January to June. 



Table 2 Available data for PG 1325+101. 



Date 


Telescope 


Filter 


Lenght 


2008/03/01 


1.5m-Loiano 


B 


5.4 h 


2008/03/13 


Im-Piszkesteto 


B 


4.8 h 


2008/02-05 


1.2 Mercator + 
Suhora & Baker 


V 


95.5 h 


2008/06/04,06 


1.6m-Moletai 


B 


3 h 


2009/01/28-31 


Im-LOAO 


B 


9 h 


2009/02/25-27 


1.2m-MONET 


B 


11.7 h 


2009/03/18 


1.2m-MONET 


B 


1.7 h 


2009/03/20-22 


2.2m-Calar Alto 


B 


11.5 h 


2009/03/25-28 
2009/03/30-31 


Im-LOAO 


B 


23.8 h 


2009/04/22 


1.2m-Loiano 


B 


2.9 h 


2009/04/29 


2.2m-Calar Alto 


B 


2.9 h 


2009/05/06 


2.2m-Calar Alto 


B 


3.9 h 


2009/05/25-27 


Im-LOAO 


B 


14 h 
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Fig. 2 0-C diagram of the main frequency (fl) of 
PG 1325+101 for the V band run. 



The pulsation frequencies of PG 1325+101 are well 
known from Silvotti et al. 2006 and Baran et al. 2010. 
We can thus fix the frequencies in order to calculate 
the mean phase of the whole data set using the soft- 
ware Period04 (Lenz and Breger 2005): this represents 
the Calculated value in the 0-C diagram. We repeat 
the procedure also for the single runs, measuring the 
Observed values. The phase differences are then turned 
into time-lag and plotted against the observations pe- 
riods. At this stage of the analysis we prefer to per- 
form the calculations using only the main pulsation 
frequency (at 7255.5 fiHz, with an amplitude of about 
27 mma), in order to avoid confusion on the sinusoidal 
fit and on the estimate of the phase with P04 for the 
runs with poor frequency resolution. 

We have calculated the 0-C diagram mantaining 
separated the two years, since we don't have enough 
data to link them. Furthermore we have the possibil- 
ity to use a long run (3 months) performed in visual 
band for another project on the same star (Baran et al. 
2010). From Figure [5] we see that the 0-C diagram of 
the V band monitoring shows a good phase coherence 
over this time scale. The uncertainies are the formal 
errors of the sinusoidal fit, according to Montgomery 
and O'Donoghue 1999. 

More difficulties are noticed when we calculate the 
0-C for the B-band data, probably related to a poor 
coverage. Figure [3] shows preliminary results for the 
three B-band runs of 2008. Considering that the second 
point has a very large error bar, the phase appears 
still rather stable. However, things are more complicate 
when we join B and V data and we are still working to 
understand what is the reason of these difficulties. 

Similar problems are found also with the data of 
2009 and preliminary results are under analysis. The 
present coverage in B band is not yet enough to ensure 
a unique solution, even over a single season. New data 
will help to converge towards a unique solution joining 
together the data of the different seasons with an im- 
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Fig. 3 0-C diagram of the main frequency (fl) of 
PG 1325+101 for season 2008 (without V band run). 

proved measurement of the pulsation periods. 
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